Ionization rates in a Bose-Einstein condensate of metastable Helium 
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We have studied ionizing collisions in a BEC of He*. Measurements of the ion production rate 
combined with measurements of the density and number of atoms for the same sample allow us 
to estimate both the 2 and 3-body contributions to this rate. A comparison with the decay of 
the number of condensed atoms in our magnetic trap, in the presence of an rf-shield, indicates 
that ionizing collisions are largely or wholly responsible for the loss. Quantum depletion makes a 
substantial correction to the 3-body rate constant. 



The observation of Bose-Einstein condensation of 
metastable helium (He in the 2 3 S\ state, denoted He*) 
jl], P constituted a pleasant surprise for experimentalists 
although the possibility had been predicted theoretically 
||. Success hinged, among other things, on a strong 
suppression of Penning ionization in the spin-polarized, 
magnetically trapped gas. Too high a rate of ionization 
would have prevented the accumulation of a density of 
atoms sufficiently high to achieve evaporative cooling. 
The ionization rate is not completely suppressed how- 
ever, and when the atomic density gets high enough, a 
cold, magnetically trapped sample of He* does produce 
a detectable flux of ions. As shown in Q, this signal can 
even be used as a signature of BEC. The observation of 
ions from the condensate opens the possibility of mon- 
itoring in real time the growth kinetics of a condensate 
||. This is an exciting prospect, but in order to quan- 
titatively interpret the ion rate, it is necessary to know 
the relative contributions of 2 and 3-body collisions. 

In this paper we use the unique features of metastable 
atoms to detect, in a single realization of a BEC, the ion- 
ization rate, the density and the number of atoms. This 
allows us to extract 2 and 3-body rate constants without 
relying on fits to non-exponential decay of the number of 
atoms, which require good experimental reproducibility 
fj. |, gland are notoriously difficult to interpret quantita- 
tively] a. After estimating the ionization rate constants, 
a comparison with the observed decay of the number of 
atoms reveals no evidence for collisional avalanche pro- 
cesses. Thus, by contrast with 87 Rb |§, He* seems to be 
a good candidate for studying "hydrodynamic" regime 
as well as the effects of quantum depletion. Indeed in 
our analysis of the 3-body ionization process, quantum 
depletion makes a substantial correction 0]. 

Much theoretical {§ |(J and experimental @, §, 0, §J 
work has already been devoted to estimating inelastic 
decay rates in He*. The dominant 2-body decay mech- 
anisms, which we will refer to collectively as Penning 
ionization, 



are known to be suppressed by at least 3 orders of mag- 
nitude in a spin-polarized sample, but the total rate con- 
stant has not yet been measured. The theoretical esti- 
mate of the rate at 1 is - 2 x 10~ 14 cm 3 s" 1 {§ |o|. 
The 3-body reaction, 

He* + He* + He* -> He* 2 +He*(~ lmK) 
^ He+ + He{lS) - 



( 2 ) 



proceeds via 3-body recombination followed by autoion- 
ization of the excited molecule. The rate has been 
estimated theoretically [[l3| to have a value of order 
1CP 26 cm 6 s _1 . Both reactions yield one positive ion 
which can easily be detected in our apparatus. 

We define collision rate constants according to the den- 
sity loss in a thermal cloud Q: ^ = — ^ — [3 n 2 — L n 3 
with n the local density, r the (vacuum limited) lifetime 
of the sample, and (3 and L the 2-body and 3-body ioniz- 
ing rate constants defined for a thermal cloud. We have 
assumed here that there are no other loss processes. One 
can calculate an expected ionization rate per trapped 
atom (r): 



r 



Ion rate 
Nn 



1 



+ - k 2 (3 n + — k 3 L n 2 , (3) 



7 



63 



He* + He* 



He+ 
He+ 



He{lS) + e- 

e~ 



(1) 



for a pure BEC in the Thomas- Fermi regime with a num- 
ber of atoms No, and a peak density n . The numerical 
factors come from the integration over the parabolic spa- 
tial profile and the fact that although 2 or 3 atoms are 
lost in each type of collision, only 1 ion is produced. The 
effective lifetime r' > r is due to ionizing collisions with 
the background gas. The factors Ki take into account the 
fact that the 2 and 3-particle local correlation functions 
are smaller than those of a thermal cloud. For a dilute 
BEC k 2 = 1/2! and k 3 = 1/3! fj §. Because the He* 
scattering length is so large, quantum depletion lead to 
significant corrections || to the k's as we discuss below. 

Much of our setup has been described previously |l], 
[l5| , |l6|| . Briefly, we trap up to 2 x 10 s atoms at 1 mK in 
a Ioffe-Pritchard trap with a lifetime (t) of 90 s. We use 
a "cloverleaf configuration [jl7j with a bias field Bo = 
150 mG. The axial and radial oscillation frequencies in 
the harmonic trapping potential are i/u — 47 ± 3 Hz and 
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V!_ = 1800 ± 50 Hz respectively (uj/2tt = (Vp^) 1 / 3 = 
534 Hz). A crucial feature of our set up is the detection 
scheme, based on a 2 stage, single anode microchannel 
plate detector (MCP) placed below the trapping region. 
Two grids above the MCP allow us either to repel positive 
ions and detect only the He* atoms, or to attract and 
detect positive ions produced in the trapped cloud. 

To detect the ion flux, the MCP is used in counting 
mode: the anode pulses from each ion are amplified, dis- 
criminated with a 600 ns deadtime and processed by a 
counter which records the time delay between succes- 
sive events. Typical count rates are between 10 2 and 
10 4 s _1 . We have checked that the correlation function 
of the count rate is flat, indicating that there is no double 
counting nor any significant time correlation in the ion 
production. The dark count rate is of order 1 s _1 . By 
changing the sign of the grid voltage, we have checked 
that while counting ions, the neutral He* detection rate 
is negligible compared to the ion rate (less than 5 %) even 
when the radio frequency (rf) shield is on. We estimate 
the ion detection efficiency by assuming that only ions 
which hit the open channels of the MCP (60% of the to- 
tal area) are detected (with a 100% quantum efficiency). 
We then multiply by the transmission of the two grids 
(0.84) 2 . Based on Refs. @, ||, we assume this (0.42) is 
an upper limit on our detection efficiency. 

To find the values of No and no corresponding to the 
measured ion rate, we use the MCP to observe the time- 
of-flight signal (TOF) of the He* atoms released from the 
rapidly switched off trap. The instantaneous count rate 
can be as high as 10 6 s _1 , and the MCP saturates when 
used in counting mode. To avoid this problem, we lower 
the MCP gain, and record the TOF signal in analog mode 
with a time constant of 400 microseconds. Several tests 
were performed to verify the linearity of the detector. 

In a typical run, forced evaporative cooling takes place 
for 40 s, down to an rf-knife frequency of 500 KHz, about 
50 kHz above the minimum of the trapping potential. 
Near the end of the ramp, the ion rate increases sharply, 
signaling the appearance of a BEC (Fig. 4 in Q|). After 
reaching the final value, the rf-knife is held on at that 
frequency (rf-shield). The above sequence results in a 
quasi pure BEC for delay times up to 15 s (see Fig. ||). 
By quasi-pure we mean that we see no evidence of any 
thermal wings in signals such as shown in the inset of 
Fig. |l|. From tests of our fitting procedure, we estimate 
that the smallest thermal fraction we can distinguish is 
about 20 %, corresponding to a temperature on the order 
of the chemical potential. Runs in which thermal wings 
were visible were discarded. 

To acquire the TOF signals corresponding to a given 
ion rate, we turn off the rf-shield, wait 50 ms, and then 
turn off the magnetic trap, and switch the MCP to analog 
mode. To be sure that the rf has no influence on the ion 
rate, we use only the number of ions observed during the 
50 ms delay to get the rate. We fit the TOF signals to an 
inverted parabola squared as expected for a pure BEC in 
the Thomas-Fermi regime, and for a TOF width (~ 5 ms) 
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FIG. 1: Chemical potential versus number of detected atoms 
to the power | and its linear fit. Data are for quasi-pure BEC. 
Inset shows a typical TOF signal and its inverted parabola 
squared fit. 



narrow compared to the mean arrival time (100 ms) [|lj. 
Under these assumptions, the chemical potential u de- 
pends only on the TOF width, the atomic mass and the 
acceleration of gravity |po|| , and thus can be measured 

quite accurately. Figure [l] shows that a varies as N^ 5 as 
expected, over almost 2 decades in atom number. Resid- 
uals from the linear fit do not show any systematic varia- 
tion which is a good indication of the detection linearity. 
A fit on a log-log plot gives a slope of 0.39. 

To determine the collision rate constants (3 and L, we 
need an absolute calibration of the number of atoms and 
the density. As discussed in Ref. M, all the atoms are 
not detected, and the direct calibration has a 50 % un- 
certainty which is responsible for the large uncertainty 
in the scattering length a. In fact the measurement of 
the chemical potential gives an accurate value for the 
product no a = [im J Anh 2 , and with the value of To gives 
the product N a = (1/15) {h/mlu) 1 ' 2 {2^/Tiidf/ 2 as well. 
Therefore, in the hopes that the He* scattering length 
will be measured more accurately in the future, we shall 
express Nq and uq in terms of a. In this paper, unless 
stated otherwise, we suppose that a = 20 nm, and in our 
conclusions we shall discuss how our results depend on a. 

Figure || shows the ion rate per atom T versus the peak 
density. The densest sample corresponds to iV = 2 x 10 5 
atoms and n = 2.5 x 10 13 cm -3 . The corresponding 
Thomas-Fermi radii are rj_ — 5 //m and rn ~ 200 /^m. 
The vertical intercept in Fig. ^| corresponds to ioniz- 
ing collisions with the background gas (1/V). We have 
independently estimated this rate using trapped thermal 
clouds at 1 mK and 5 /iK, and found 1/V < 5 x 10~ 3 s _1 . 
This value is negligible at the scale of the figure. 

The curvature seen in Fig. || shows that 3-body ioniz- 
ing collisions are significant. Before fitting the data to get 
(3 and L, we must take into account several effects. First, 
for 3-body collisions, quantum depletion is important. 
For T = 0, reference || gives a multiplicative correction 
|l| to the factor k 3 of ((1 + e) = (1 + A-p| x Vn~o~a*)), 
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FIG. 2: Ion rate per trapped atom versus peak density for 350 
different quasi- pure BEC's. Atom number and density are 
deduced from /i , u and a (here 20 nm). Data were taken for 
2 different bias fields corresponding to v± = 1800 Hz (crosses) 
and f± = 1200 Hz (circles). The dashed line corresponds to 
the best fit involving only 2-body collisions. The solid line is 
a fit to 2 and 3-body processes. 



where A ~ 0.84 and comes from an integration over the 
spatial profile using a local density approximation. At 
our highest density e ~ 0.5. Two-body collisions are 
subject to an analogous correction but approximately 3 
times smaller. The fits in Fig. || include the density 

dependence of ^2,3, associated with quantum depletion. 

3/2 

The n dependence introduced for 2-body collisions is 
far too small to explain the curvature in the data. The 
density dependence of «2,3 does not significantly improve 
the quality of the fit, but it significantly affects the value 
of the fitted value of L (reduction of 40%). 

In addition, the fact that the sample probably contains 
a small thermal component means that collisions between 
the condensed and the thermal parts must be taken into 
account ||, Assuming a 10 % thermal population 
(j^f — 1-1), we find K3 = g(l+e+e'), with an additional 
correction e' ~ 0.35 for the densest sample p2| . 

Taking into account all these corrections, the fitted 
values of the collision rate constants are: P20 — 
2.9(±2.0) x 10- 14 cm 3 sec" 1 and L 20 = 8.5(±5.3) x 
10~ 27 cm 6 sec -1 , where the subscripts refer to the as- 
sumed value of a. These values are in good agreement 
with the theoretical estimates. The error bars are esti- 
mated as follows. We fix either (3 or L and use the other 
as a fit parameter. We repeat this procedure for differ- 
ent values of the fixed parameter and take the range over 
which we can get a converging and physically reasonable 
fit (i.e. no negative rate constants) as the uncertainty 
in the fixed parameter. These error bars are highly cor- 
related since if is increased, L must be decreased and 
vice-versa. The error bars do not include the uncertainty 
in the absolute ion detection efficiency (see below). 

Until now we have assumed a = 20 nm, but current ex- 
periments give a range from 8 nm to 30 nm [jj], ^) . Using 
Eq. 3 and our parameterization of no and Nq in terms 
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FIG. 3: Fraction of remaining atoms measured by TOF as a 
function of time. The rf-shield is on and the cloud remains a 
quasi pure condensate during the decay. The lines correspond 
to the predicted atom decay according to Eq. 3 with the 
fitted value of the 2- and 3-body rate constants for a — 10 nm 
(dashed line), a = 20 nm (solid line) and a = 30 nm (dotted 
line). The case of a = 10 nm is not necessarily excluded 
because other, non-ionizing losses could be present. 



of a, one can see that, in the absence of quantum deple- 
tion, the values of [3 and L extracted from our analysis 
would be proportional to a 2 and a 3 respectively. Taking 
quantum depletion into account, no simple analytical de- 
pendence exists, but one can numerically evaluate (3 and 
L vs. a and fit the results to expansions with leading 
terms in a 2 and a 3 respectively. The effect of quantum 
depletion is negligible for (3 {(3 a ~ P2o{^) 2 )- For L, we 
find L a w L 20 (|j) 3 [l - 0.23^-P] with a in nm. 

To test the consistency of our measurements, we have 
studied the decay of the number of atoms in the BEC 
(Fig. ||). To acquire these data, we held the BEC in the 
trap in the presence of the rf-shield for varying times. 
This study involves multiple realizations of a BEC, which 
typically exhibit large fluctuations in the initial atom 
number. We have been able to reduce this noise by us- 
ing the ion signal to select only data corresponding to the 
same ion rate 500 ms after the end of the ramp. This time 
corresponds to t = in the figure. We have also plotted 
the predicted decay curve (solid line) corresponding to 
ionization only. This curve results from a numerical in- 
tegration of the atom loss due to ionization processes, 
calculated from the fitted values /3 2 o an d £20- The fact 
that the error bars on (3 and L are correlated leads to 
a small uncertainty on the solid curve that happens to 
be of the same order of magnitude as the typical error 
bars on the data. The observed decay agrees fairly well 
with the solid curve, and ionization apparently accounts 
for most of the loss. If the ion detection efficiency were 
actually lower than we assume, the predicted decay due 
to ionization would be faster than the observed decay, an 
unphysical situation. From this, we conclude that our 
estimate of the detection efficiency is accurate and does 
not lead to an additional uncertainty in (3 and L. 
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We have also plotted the curves obtained from the 
same analysis but with scattering lengths of 10 and 
30 nm. The curve corresponding a = 30 nm lies below 
the data points. Based on our analysis, this means that 
a = 30 nm is excluded. A scattering length of a = 25 nm 
is the largest one consistent with our data. In contrast, 
the decay predicted for an analysis with a = 10 nm is 
slower than the observed decay of the number of con- 
densed atoms. This would mean that there are additional 
non-ionizing losses (contributing up to half of the total 
loss), and/or that we have overestimated the ion detec- 
tion efficiency by a factor as large as 2. In the latter case, 
the rate constants [3 and L should be multiplied by the 
same factor. This results for a = 10 nm in a supplemen- 
tary systematic uncertainty on (3 and L of a factor as 
large as 2. 

Even though the peak densities of our BEC are small 
compared to those in alkalis, the elastic collision rate 
is high because of the large scattering length, and one 
must consider the possibility of collisional avalanches. 
For a = 20 nm our densest cloud has a mean free path of 
(87ra 2 n) _1 ~ 7fim and using the definition of || the colli- 
sional opacity is w 0.8. With Rb atoms this would result 
in avalanche processes increasing by a large amount the 
atom loss § . In our case we have to consider secondary 
collisions leading to both ion production and atom loss. 
However, no product of an ionizing collision can pro- 
duce many secondary ionizing collisions, since the cor- 
responding mean free paths are at least two orders of 



magnitude larger than rii . Hence secondary ionization is 
unimportant and Eq. [| correctly describes the ionizing 
processes. This conclusion is supported by our observa- 
tion that no correlation exists in the time distribution of 
detected ions. 

The good agreement between the data and the curve in 
Fig. 3 indicates that losses due to non-ionizing collisional 
avalanches are not taking place either. This is in agree- 
ment with data on elastic collisions with He + , He^ and 
He(lS'), which have small cross sections |2j|. Collisions 
with hot He* atoms from the reaction of Eq. || arc more 
likely to play a role, but due to the higher velocity, the 
elastic cross section for these atoms is smaller. This is in 
contrast to what happens with Rb atoms || where the 
total cross section is enhanced due to d-wave scattering 
resonance. 

The theoretical analysis shows that quantum deple- 
tion strongly affects the measured 3-body rate constant. 
One way to experimentally demonstrate this effect would 
be to do similar measurements with thermal clouds, and 
compare them with the results reported here. Absolute 
calibration of ion and atom detection efficiency should 
play no role in this comparison, if one could prove that 
they are the same for both situations. 
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